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Effects of order parameter collective modes on electromag- 
netic response are studied for a clean spin-triplet supercon- 
ductor with k x ± iky orbital symmetry, which has been pro- 
posed as a candidate pairing symmetry for Sr2Ru04. It is 
shown that the k x ± ik y superconductor has characteristic 
massive collective modes analogous to the clapping mode in 
the A-phase of superfluid 3 He. We discuss the contribution 
from the collective modes to ultrasound attenuation and elec- 
tromagnetic absorption. We show that in the electromagnetic 
absorption spectrum the clapping mode gives rise to a reso- 
nance peak well below the pair breaking frequency, while the 
ultrasound attenuation is hardly influenced by the collective 
excitations. 
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Superconducting and superfluid properties of uncon- 
ventional Cooper pairing states continue to be of interest 
in the research of strongly correlated fermion systems 
such as liquid 3 He, heavy fermion compounds, high-T c 
superconductors and, the more recently discovered super- 
conductor Sr2Ru04tl The spin-triplet k x ±ik y -wave state 
is an example of the unconventional pairing states and 
is known to describe the superfluidity of 3 He- A. Two- 
dimensional (2D) version of this state has recently at- 
tracted much attention as a candidate pairing state for 
Si'2Ru04. 

The possibility of the spin-triplet p-wave pairing in 
Sr 2 Ru04 was first suggested by Rice and SigristB from 
the analogy to 3 He. There now exist considerable exper- 
imental evidences supporting the unconventional spin- 
triplet superconductivity in SrjRu04. The absence of 
the Hebel-Slichter peak in l/TiB and substantial reduc- 
tion in T c by non-magnetic impuritiesu indicate that 
Si'2Ru04 is at least not a conventional s-wave supercon- 
ductor. The fj.SK experimental suggest a pairing state 
with broken time reversal symmetry. More strong evi- 
dence of the triplet pairing has been found by 17 NMR 
measurements^] in which the temperature-independent 
Knight shift is observed for magnetic fields parallel to 
the RuC>2 plane. Possible pairing states in Sr2Ru04 have 
been classified from a group theoretical point of view.EI 
It was shown that the quasi-2D electronic structure of 
Sr 2 Ru04 leads to five possible p-wave states stabilized 
by the absence of gap nodes. Among these states, the 
pairing state compatible with all the above experiments 
is the k x ± iky state which is the 2D analog of 3 He-A. 

Study of the electromagnetic (EM) response provides 
valuable information on the properties of unconventional 



superconductors. Dynamical properties are, in particu- 
lar, intriguing because they depend not only on the equi- 
librium gap structure but also, on the collective excita- 
tions of the order parameter .Em As is known in the study 
of superfluid 3 He, internal degrees of freedom of the or- 
der parameter give rise to the order parameter collective 
modes (OPCM's) specific to a given pairing symmetry. 
In this paper, we discuss, within the collisionless regime, 
how the OPCM's in the k x ± ik y superconductor con- 
tribute to the EM response. 

The OPCM's in the, 2D k x ± ik v state has been dis- 
cussed by TewordtBII2l As in the case of 3 He- A, there 
exists the clapping mode, which is a characteristic mas- 
sive collective mode resulting from the orbital degrees of 
freedom of the k x ± ik y -wave order parameter. Tewordt 
showed that the coupling of this mode to external fields 
vanishes in the long wavelength limit (q — > 0)t3 and so 
he did not attempt to discuss the observability of the 
clapping mode. The effect of a finite wave vector q is, 
however, important even in type II superconductors when 
we consider the role of OPCM's in the EM response. We 
in fact show in this paper, by calculating a g-dependent 
dynamical conductivity, that the clapping mode leads to 
an EM absorption peak. We also discuss the ultrasound 
attenuation. We show that, although ultrasound experi- 
ments have played a key role in the study of the OPCM's 
in 3 He, the ultrasound cannot be a good probe of the 
OPCM's in metals such as Sr2Ru04. This is because the 
sound velocity in ordinary metals is much smaller than 
the Fermi velocity v-p, in contrast to 3 He. 

The 2D k x ± ik y state is defined on a cylindrical Fermi 
surface and the orbital structure of the order parameter 
is expressed in terms of two basis functions k x — cos 9 
and k y — sin 6 specifying the direction of the Fermi mo- 
mentum kp. The matrix order parameter is given by 



A k = A(k x ± ik y )u z ia y — Ae a 



(1) 



where <jj's are Pauli matrices. The gap (AkA^) 1 / 2 = A 
of the two unitary states is independent of k and the 
Bogoliubov-quasiparticle energy = + AkA^) 1 / 2 is 
the same as that of the s-wave state. 

The collective excitations of the order parameter are 
described as oscillations of order parameter fluctuations 
from the equilibrium state. In the k x ± ik y state, one of 
the two degenerate states is realized in equilibrium. We 
choose the k x + ik y state as the equilibrium state and 
consider an order parameter fluctuation of a plane wave 
form (5Ak(q, cj)e^ q ' r ~ wt ) . Such fluctuation can be excited 
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by space-time dependent EM fields, a scalar potential 
0(q, Lo)e^ r -^ and a vector potential A(q, w)e l(q r " wt) . 
The order parameter fluctuation is expressed using the 
two bases e ±jfl 



<JA k (q,u;) - [D+(q, ^)e iy + £>_(q, wje""]^*?-,,, (2) 

where D + and _D_ are variables to be determined by the 
self-consistency equation. Since D + and Z?_ are complex 
variables, there are four degrees of freedom of the order 
parameter fluctuations in the representation of (|2|). 

A well-established way to study the dynamical prop- 
erties including the fluctuations of the order param- 
eter is to introduce matrix distribution functions in 
spin space, n k (q, <jS) and /k(q, u>) with matrix ele- 



ments (n k )a/: 
Jdte 



Jdt, 



\ a k_/3 a k; + a/ 



and (/ k ) 



a/3 



{a^ k0 a 



k J (k± = k± q/2). The deviation 
<5/k of f k from its equilibrium value determines the order 
parameter fluctuation via the self-consistency equation: 



6A k (q, u) 



^\> k , k '<5/ k '(q,w), 



(3) 



(£ k+ -E k _)(e k _ -e k , 

oj* - (E k+ - E k _)* 



(7) 



Here we have dropped the arguments q and to for brevity 
(6A 1 * = 5A^(-q,-cj)). In Eqs. (|) and (§), small correc- 
tions of order q/kp have been neglected and accordingly 
A k± has been put as A k .liil 

Now we discuss the OPCM's in the k x + ik y state. It 
is convenient to introduce the following new variables: 



l[D + (c L ,oj)±D* + (- q ,-u;)] = S^ 



-''"•±D*(-q,-w)e 2ifl «] = ' D ' 



D" 



(8) 
(9) 



where e l9q = q x +iq y . Using Eq. (|J), the self-consistency 
equation @ can be written in terms of these variables as 



Mi 



M 2 



D" 





2zA{^ r (o;5eo + ?/£ei) sin 2d) k 

( 2A(T{ude + v5ei)) k 
\ 2A(F(ujSe + rjSex) cos 2i?) k 



(10) 



(11) 



where u kik < = —2g\k ■ k' is the pairing interaction. The 
deviation 5n k is related to the charge current density as 



Ti 2 

J(q,w) = ev F VkTr<5n k (q,w) A(q,w), (4) 

' — ' mc 



where n is the number density. (Note that for the spin- 
independent perturbations considered here, Sn k is pro- 
portional to the unit matrix and thus the trace in Eq. 
(f|) gives only a factor 2.) The distribution functions have 
been intensively studied in the context of ultrasound at- 
tenuation in superfluid 3 He (see, for a review, Ref. 12). 
Assuming the particle-hole symmetry, we can write the 
distribution functions integrated over the energy variable 
6c as 



dtf k (5n k = (Se + 8ex)r)/ (u - i]) 



- 2^A k Aj i (w ( 5eo + v5ei)/(uJ ~ v) 
+ ^(cu + i 1 )(6A k At-A k 6At), 

d&6f k = r d&^6A k 

J-e, E k 



(5) 



T 



rl 



^( W 2 -2A k A k -ry 2 )(5A k 
-A k 5A k A k - (a;5eo+77<5ei)Akl, (6) 



where SeQ = e(f>(q,uj) and Se\ — — (e/mc)kp • A(q, uj) 
represent the perturbation energies, i] = vp • q, k = 
— (1/2) tanh(i?k/2T), e c is the usual cutoff energy and 



T = 



'-oo ^E k E k _ 



(E k+ + £ k J(9 k _ 



- (E k+ + E k 



where $ is the angle between k and q, (• • -) k denotes the 
angle average over the Fermi surface and Mi j2 are 2x2 
matrices with the following elements: 



V])k, 



{M 1 ) ll = {T(u 2 -AA 2 ~ V 2 )) kl 
(Mi) 22 = (T\u? - 2A 2 (1 + cos4tf) 
(Mi)ia = (Mi) a i - {F{lj 2 - 4A 2 - 7 ? 2 )cos2z9) l 
(M 2 ) n - (T{uj 2 - v 2 )) k , 
(M 2 ) 22 = {T[uj 2 - 2A 2 (1 - cos4tf) - tj 2 }^ 
(M 2 )ia = (M 2 ) 21 = {T(uj 2 -77 2 )cos2^)k. 



(12) 



Note that the variables D' and D" are decoupled from 
each other. 

The OPCM's correspond to the eigen modes in the 
absence of external fields. The eigen frequencies u> are 
therefore determined by detMi. 2 = 0. In the long wave- 
length limit, q — 0, it is easy to obtain the eigen frequen- 
cies. In this limit, since the off-diagonal elements of M\p, 
may be neglected, all the modes decouple. The results for 
q — > are summarized in Table [jj. The D^_-mode has the 
eigen frequency just at the pair breaking edge, u = 2A. 
The variable Z?" represents a mode corresponding to the 
phase fluctuation of the order parameter, namely, the 
so-called Anderson-Bogoliubov mode. This is a gapless 
mode with sound- like dispersion ui = {v-p/\2)q but is re- 
placed by a plasmon when the Coulomb interaction be- 
tween electrons is taken into account .0 The remaining two 
modes, £>_-modes, have the same character as the clap- 
ping mode in superfluid 3 He-A.E3ii3 The eigen frequency 
of the clapping modes in the 2D system is u> = \/2A and 
is below the pair breaking energy 2 A. At finite q, the 
dispersion relation of the clapping modes up to the order 
q 2 is given by 



\(q) = 2A 2 + hv F q) 



(13) 



2 



The coupling of the OPCM's to the external fields is 
determined by the right-hand side of Eqs. (|Io| ) and (11). 
We see that the -D±-modes and the D^-modes couple 
to the transverse field and the longitudinal field, respec- 
tively. 

The presence of the coupling between the D^.-modes 
and the longitudinal field means that these modes can 
be excited by a longitudinal phonon field. Since the 
collective excitations couple to the density fluctuation 
Sp(q,u>) = ^ k Tr<5rik via the last term in Eq. (||), there 
is a possibility that the clapping mode (D"-mode) is 
observed by ultrasound measurements. We have esti- 
mated the ultrasound attenuation coefficient by assuming 
typical conditions u> <C v^q <C A for ultrasound mea- 
surements in ordinary metals and also in S*«Ru04. (In 
Sr2Ru04, the sound velocity is ~ 10 5 cm/alJ, the Fermi 
velocity is ~ 10 7 cm/sB'El and A - 1 K.) The attenu- 
ation coefficient from quasiparticle excitations coincides 
with that for the s-wave superconductor because of the 
isotropic gap. The contribution from the collective ex- 
citations can be obtained using Eqs. (§) and ([ll]). Tak- 
ing into account the above conditions and considering 
the limit uj — > 0, one can find, after straightforward al- 
gebra, that the collective excitations give only a small 
correction of order (uj/v-pq) 2 ~ 1CP 4 to the quasiparti- 
cle contribution. Consequently, the ratio of the super- 
conducting to normal attenuation coefficient in the low- 
frequency limit ui — > is given by the well-known BCS 
result a s /a n — 2/(e A / T + 1). Thus sound waves are not 
suitable probe for the OPCM. It is to be noted that, in 
contrast to ordinary metals, ultrasound propagation in 
superfluid 3 He is in the high frequency regime such that 
ui ~ A ^> vpq where the OPCM plays an important role. 
This is a reason why sound waves can be a good probe 
of the OPCM's in superfluid 3 He. 

Now we consider the transverse response of a k x + ik v 
superconductor. The transverse response is conveniently 
described by the complex conductivity at (q, ui) for the 
transverse electric field E t (<\,uj) = (iuj/c)A t (q, uj). The 
real part o\ of the complex conductivity at = o~\ + «cr 2 
determines the absorption of EM waves. We can readily 
obtain at by calculating the transverse current density 
using the above formulation. The conductivity consists 
of three characteristic parts, i.e., a t = of p + cr™ 1 + af la ; 
the first and the second terms are contributions from 
quasiparticle excitations and from collective modes, re- 
spectively, and the last term of la = ine 2 /mui arises from 
the diamagnetic current (— ne 2 /mc)At. The quasiparti- 
cle contribution a^ p is the same as the BCS result for the 
s-wave superconductor: 



cr t cm (q,w) = -2iev F A^2k t r)Tsm2tiD'_/E, 

k 



T qp 



(q,w) 



2ine 2 k 2r q 2 



x ui 2 — r\ 2 



;i - 2A 2 ^)) k , 



(14) 



where k t = sin $ is the transverse component of the unit 
vector k. The collective-mode contribution cr™ 1 comes 
only from the D'_ -clapping mode and is obtained as 



4me 2 A 2 (T(uj 2 - 4 A 2 - r] 2 )) k (Tr]kt sin2z9) 2 



muj 



detMi 



(15) 



where the last line is derived using D'_ determined by 
Eq. |[C]). Note that Eq. (|l|) includes detMi in the de- 
nominator. This implies that a resonance EM absorption 
due to the clapping mode occurs at the eigen frequency. 
To see this explicitly, it is useful to estimate the real 
part <7i by assuming the long wavelength condition such 
that v-pq <C A. This assumption is justified for type II 
superconductors, since the penetration depth A is large 
compared with the coherence length £o and the important 
values of q in the Fourier integral of the fields penetrating 
into the superconductor are < A . In what follows, we 
discuss the w-dependence of a\ in this case. 

Since we are interested in high frequencies ui ~ A, we 
may expand the conductivity in terms of q. Expanding 
Eqs. ([Ij) and ( jl5| ) up to the order q 2 and adding the 
results, we obtain 



0t(q,a;) = 



2 2 2 

me vf,q 
Amu) 3 



l-T 



4A' 



(16) 



where To — A 2 T(q = Q,lu). Equation Jig 1 ) clearly shows 
that at has a pole at the eigen frequency of the clapping 
mode, yielding a delta function peak in the real part <7i . 
The diamagnetic term af la in Eq. (|l6| ) is pure imaginary 
except at uj — 0, so that only the first term contributes 
to a\ at the high frequencies of interest. 

We briefly discuss a\ in the low frequency region 
uj ^ A where the OPCM does not come into play. 
There are two contributions to a\ at the low frequen- 
cies. One is a delta function at ui = 0, associated with 
superfluid motion excited by EM fields. The other is 
caused by thermally excited quasiparticles. The energy 
conservation u> = |£' k+ — Ek_ \ for the corresponding 
quasiparticle-scattering process requires that this contri- 
bution is restricted in the region uj < v-pq <C A. We note 
that in the limit vpq/A — > 0, the conductivity sum rule 
f divai(q,uj) = une /2m is satisfied only by the two 
contributions. As can be shown from Eq. ( |l4|) , in the 
limit v-Fq/A — > 0, the ^-integration yields irri^e 2 /2m for 
the S(uj) term and irn n e 2 /2m for the contribution from 
thermally excited quasiparticles, where n s is the super- 
fluid density and n n — n — n s is the normal- fluid density. 
This explains why di at w ~ A, given by Eq. (|l6|), is 
proportional to q 2 . 

Equation ( |l6| ) depends on temperature T through To 
and A. At sufficiently low temperatures, the hyperbolic 
tangent function, tanh(i?k/2T), in To may be replaced 
by unity; then Eq. ( |l6|) is reduced to the result for T = 0. 
Figure |l| shows the w-dependence of the real part a\ at 
T = obtained from Eq. ([!(]). The structure above uj = 
2 A is caused by pair-breaking processes. Below the pair- 
breaking edge, uj < 2A, there is a delta function peak at 
uj c \ = V2A, associated with the absorption of the clap- 
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ping mode. This collective-mode contribution is explic- 
itly given by erj m (chw) = (ir 2 ne 2 Vpq 2 /64mA 2 )S(uj - w cl ). 

Hirschfeld et alB have studied the power absorption 
from the OPCM in the Balian-Wcrthamer state (in a 3D 
system). The power absorption P(uj) is given as the dissi- 
pation of the field energy which appears as Joule heat. To 
consider -P(w), we need to take into account the presence 
of the vacuum-metal interface. They estimated P(io) by 
assuming the specular surface boundary condition and by 
neglecting the surface scattering effect causing the pair 
breaking near the surface. In the 2D k x ± ik y state con- 
sidered here, the same calculation is repeated by using 
Eq. Jl6| ) and considering the EM wave injected along the 
2D plane. The resulting collective-mode contribution to 
P(lu) has the following w-dependence: P cm (w) = be- 
low u> = w c i(0). Above the threshold, u> > w c i(0), P cm (u) 
increases rapidly with to and then has a peak structure 
with a finite width ~ (£ /Al)A (A l = (mc 2 /47me 2 ) 1 / 2 , 
the London penetration depth); this (^-dependence is de- 
scried by P cm (w) ~ [w 2 -w 2 (0)] 1 / 2 /[o; 2 -w 2 1 (0)+c 2 A^ 2 ] 2 , 
where c\ — vf/V% is the velocity determining the dis- 
persion of the clapping mode. In the power absorption 
spectrum, the OPCM gives rise to not a delta function 
peak but a broadened peak. This is a consequence of 
taking into account the dispersion of io c \ in the denomi- 
nator of Eq. (plf). In actual metals, impurity scattering 
also brings about the broadening. The power absorption 
spectrum, however, has a definite peak structure due to 
the collective excitation if £o is small enough compared 
with Xl and at least the superconductor is clean. This 
demonstrates that the clapping mode in the k x ±ik y -w&ve 
type II superconductor can be observed by EM absorp- 
tion measurements. 
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FIG. 1. Real part eri(q, ui) of the complex conductivity at 
T = as a function of reduced frequency lo) A. The ordinate 
is scaled by ne 2 Vpq 2 /mA s . A delta function peak at the eigen 
frequency of the clapping mode, lo/A — \/2, is plotted as a 
Lorentzian with a tiny width. 



TABLE I. Order parameter collective modes of the k x +ik y 
state in the long wavelength limit q — > 0. 



Variable 


Eigen frequency 


Mode 


D'+ 


2A 




D'i 





Anderson-Bogoliubov mode 


D'_, D'I 


V2A 


Clapping mode 



4 



0.05 



3 




co/A 



